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programs  generated  by  Mr.  Kerb  Klimach. 
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and  derived  the  programs  for  generating  the  magnetic  data  tapes. 
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1 

INTRODUCTION 

This  is  the  final  report  on  the  clutter  statistics  task  on 
Contract  DAAK40-77-C-0112 , Distributed  Spaceborne  Radar  Study.  On 
this  task,  Amendment  P0002  to  the  contract,  the  Environmental  Research 
Institute  of  Michigan  (ERIM)  with  the  assistance  of  Decision-Science 
Applications,  Inc.  (DSA)  as  subcontractor  obtained  the  radar  clutter 
distributions  of  L-band  data  collected  on  recent  flights  of  the  ERIM 
X-L  band  SAR  (synthetic  array  radar)  system  and  generated  clutter 
tapes  by  means  of  a radar  signal  simulation.  The  principle  tasks 
included: 

1.  SAP  image  data  digitization, 

2.  SAR  image  data  calibration, 

3.  obtaining  radar  clutter  cross-section 
distributions, 

4.  generating  clutter  models, 

5.  preparing  radar  system  models,  and 

6.  generating  clutter  tapes,  and  synthesizing 
target  signals. 

1.1  BACKGROUND  AND  OBJECTIVE 

The  objective  of  this  program  was  to  produce  clutter  distributions 
of  three  different  terrain  types  at  three  different  elevation  angles 
and  at  two  different  resolutions.  These  clutter  statistics  wojld  be 
used  as  realistic  inputs  for  radar  clutter  models. 

A meeting  was  held  at  ERIM  with  representatives  from  RADC,  DSA, 
and  ERIM  to  determine  the  values  of  the  parameters  described  above. 

As  a result  of  this  meeting,  the  L-band  horizontal  transmission- 
horizontal  receive  data  collected  by  ERIM,  of  Labrador,  Canada,  during 
February  and  March,  1977  was  selected  as  the  data  set.  Data  at  20  , 
30°,  and  60°  depression  angles  would  be  used.  An  angular  interval  of 
2°  centered  on  these  3 depression  angles  would  be  digitized.  The 
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scenes  selected  were  to  include  mountains,  ice,  and  a city.  The 
data  would  be  analyzed  at  3 meter  resolution  and  30  meter  resolution. 

The  general  approach  used  was  to  form  a radar  image  using  the 
optical  processor.  This  optical  image  would  be  converted  to  a digital 
format  using  ERIM's  image  dissector  facility.  The  dissected  and  digi- 
tized image  would  be  stored  on  CCT's  in  a format  suitable  for  statis- 
tical analysis  using  the  Univeristy  of  Michigan's  computer.  After  the 
numerical  analysis  had  been  performed,  the  numbers  stored  on  tape  would 

be  converted  to  a values. 

o 

1.2  SUMMARY  OF  REPORT 

This  report  is  divided  into  6 sections  (Volume  I)  and  3 appendices 
(Volume  II).  This  first  section,  the  Introduction,  briefly  describes 
the  objectives  of  the  program  and  the  methods  used  to  satisfy  these 
objectives. 

Section  two,  Test  Data,  of  the  report  describes  the  data  selected 
for  analysis  and  how  the  data  was  converted  from  optical  images  to 
clutter  statistics.  Included  in  section  two  are  digitization  and 
calibration  methods. 

Section  three.  Statistics  for  Terrain  Clutter,  describes  the  data 
analysis  and  discusses  the  distributions  used  to  fit  the  data.  Exam- 
ples of  the  density  distribution  and  cumulative  distributions  for  various 
terrain  are  given  in  this  section. 

Section  four,  Simulation  of  Radar  Signals  from  an  R,  R Clutter 
Map,  taken  from  a technical  memorandum  by  Philip  Tomlinson  of  DSA, 
describes  a method  of  generating  simulated  radar  signals  using  clutter 
statistics. 

Section  five.  Summary,  summarizes  the  results  and  compares  the 
results  obtained  from  this  program  to  results  published  in  the 
literature.  Section  six  contains  the  Conclusions  and  Recommendations. 
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Appendix  I consists  of  copies  of  the  computer  print-out  of  the 
statistics  generated  on  this  program.  Appendix  II  is  a copy  of  the 
computer  program  and  is  included  to  make  the  documentation  complete. 
Appendix  III  is  a copy  of  the  IPL  logs  associated  with  each  data  file 
generated  during  this  program.  Its  function  is  to  allow  the  data  user 
the  opportunity  to  retrace  the  data  back  to  its  origin,  namely  the  radar 
imagery  archived  at  ERIM.  Thirty-seven  different  data  files  were 
generated  on  this  program  and  are  stored  on  9 CCTs  in  ERIM's  radar 
image  library. 
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2 

LABRADOR  TEST  DATA 


This  section  describes  the  test  sites,  shows  examples  of  radar 
imagery  of  these  test  sites,  and  discusses  the  digitization  and  cali- 
bration procedures  used  and  the  computer  program  which  calculated  the 
statistics.  The  radar  imagery  u&ctf  in  this  program  was  collected  by 
ERIM  during  Feb  iry  and  March  of  1977  as  part  of  a program  called 
"SAR-77"  and  t msored  bv  the  Canadian  Cold  Oceans  Resources  Engi- 
neering Laboratory,  Memorial  University,  St.  John,  Newfoundland. 

During  SAR-77,  four-channel  imagery  was  obtained  from  Lake  Erie  along 
the  St.  Lawrence  River  to  the  Gulf  of  St.  Lawrence,  north  to  Goosebay 
and  Hopdale,  Labrador.  The  imagery  continues  from  Goosebay  to  St. 
John,  Newfoundland;  from  St.  John,  Newfoundland  to  Quebec  City  via  New 
Brunswick;  and  finally,  returning  to  Lake  Erie  along  the  St.  Lawrence 
River.  A map  of  this  route  is  shown  in  Figure  1.  Each  image  swath 
appears  in  four  separate  output  films  corresponding  to  the  four 
channels  of  radar  data.  These  channels  are: 

X-band  parallel  polarized, 

X-band  cross-polarized, 

L-band  parallel  polarized,  and 
L-band  cross-polarized. 

The  majority  of  the  data  was  collected  using  horizontal  transmission 
although  some  data  is  available  -sing  vertical  polarization  for 
transmission.  Although  four— channel  imagery  is  available,  only  the 
L-band  parallel  (horizontal  transmit,  horizontal  receive)  imagery  was 
used  in  this  program. 

2.1  DATA  MEASUREMENT  AND  RECORDING  METHOD 

Radar  data  is  collected  on  board  the  imaging  aircraft  and  stored 
as  Doppler  phase  histories  on  photographic  film.  This  film  is  called 
a signal  film  and  is  photographically  developed  after  the  aircraft 
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Figure  1.  X-L  Band  Radar — Complete  Route  Imaged 


lands. 


This  film  transnarency  is  used  as  the  input  to  the  optical 
processor.  The  optical  processor  converts  the  Doppler  histories  into 
a radar  image  which  normally  is  recorded  on  70mm  photographic  film 
called  the  image  film.  On  the  image  film,  intense  echo  signals  appear 
as  high  photographic  density  areas  while  weak  radar  ec  ho  components 
are  low  density  areas. 

Because  one  of  the  program  goals  was  to  generate  clutter  statis- 
tics, the  image  was  converted  into  a digital  tape  for  use  in  a digital 
computer.  This  is  accomplished  using  a device  called  an  image  dissector. 
Similar  to  the  optical  image  processor,  the  dissector  converts  data 
from  the  signal  film  to  an  optical  image.  Instead  of  storing  this  image 
on  photographic  film,  it  is  (at  this  point)  converted  into  numbers  pro- 
portional to  radar  reflectivity  (or  square  root  of  reflectivity)  for 
each  image  element  selected.  These  numbers  are  stored  on  computer 
compatible  tape  and  are  used  as  inputs  to  a digital  computer  to  gener- 
ate clutter  statistics.  This  procedure  is  shown  in  block  diagram 
form  in  Figure  2. 

2.1.1  FADAR  IMAGERY  OF  SELECTED  TEST  SITES 

The  ice  imagery  digitized  was  from  the  Labrador  Sea  and  was 
selected  because  it  contains  both  homogeneous  pack  ice  and  pack  ice  made 
nonhomogeneous  by  including  large  sections  of  refrozen  leads.  The 
refrozen  leads  are  classified  as  black  ice  because  of  their  appearance. 

A radar  image  of  the  pack  ice  is  shown  in  Figure  3.  The  areas  within 
the  narrow  strips  are  the  areas  digitized  and  used  in  the  analysis. 

As  can  be  seen  from  Figure  4,  there  are  many  frozen  mountain 
lakes  in  this  selected  scene  of  mountain  terrain.  During  the  analysis 
(section  3) , it  was  noticed  that  the  mountain  data,  at  30°  depression 
angle  had  a o about  10  dB  below  the  o at  22  and  60  . Ipon  examina 
tion  of  the  imagery,  it  was  noticed  that  most  of  the  data  at  30  was 
shadowed  by  a mountain  peak.  A new  area  at  30  depression  angle  was 
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Figure  2.  Data  Flow  Diagram 
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Figure  3.  Radar  Image  of  Pack  Ice 
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selected  and  analyzed  and  its  o value  is  more  in  agreement  with  the 

o 

values  at  22°  and  60°.  This  area  is  shown  in  Figure  5. 

The  city  of  Toronto  was  used  for  the  analysis  of  clutter  statis- 
tics of  cities.  The  entire  scene  shown  in  Figure  6 was  digitized.  Only 
the  areas  included  in  the  narrow  strips  were  used  for  analysis.  Data 
from  the  large  scene  is  available  for  future  large  area  analysis. 

The  data  from  three  "fields"  in  Ontario  were  digitized  and 
analyzed  to  provide  homogeneous  data  that  would  yield  statistics  which 
would  be  easier  to  interpret.  Also  its  averaging  characteristics 
shoul  be  better  behaved  than  some  of  the  nonhomogeneo as  clutter  data. 
These  fields  subtend  a small  angular  dimension  in  elevation  so  data 
were  only  digitized  for  one  elevation  angle.  Radar  images  of  the 
three  fields  are  shown  in  Figures  7,  8,  and  9. 

2.1.2  DIGITIZATION  OF  RADAR  IMAGES 

As  mentioned  in  section  2.1,  the  ERIM  image  dissector  facility 
[1]  was  used  to  convert  the  radar  images  of  various  terrain  types  into 
computer  compatible  tapes.  These  CCT's  were  then  used  to  compute  the 
clutter  statistics  on  the  University  of  Michigan's  Amdahl  470  computer. 
The  software  used  to  compute  the  statistics  was  developed  on  this 
program  and  is  described  in  section  2.2.2. 

The  images  used  in  this  program  were  processed  to  a resolution 
cell  of  3 meters  in  range  and  4.5  meters  in  cross  range  or  azimuth 
but  was  sampled  by  the  image  dissector  using  a 1.5  meter  by  1.5  meter 
pixel  size.  This  was  accomplished  by  adjusting  the  optical  magnifica- 
tion between  the  optical  processor  and  the  image  dissector  so  that  the 
radar  azimuth  resolution  elt’  ?as  twice  the  size  of  the  dissector 

azimuth  resolution  element  10. 

To  study  the  effect  o terent  resolution  on  the  clutter  statis- 

tics, different  quantities  of  data  were  averaged  during  the  statistical 
computation  process.  All  the  initial  image  data  was  processed  to  3 m by 

1.  D.  A.  Ausherman,  "Radar  Data  Processing  and  Exploitation  Facility, 
Record  from  the  IEEE  International  Radar  Conference,  1975. 
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Figure  5.  Radar  Image  Mountainous  Terrain, 
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Figure  6.  Radar  Image  of  Toronto,  Ontari 


igure  9.  Radar  Image  of  Field 
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Figure  10.  Radar  Resolution  Cell 
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4.5  m resolution.  The  resolution  of  the  clutter  patches  was  then 
varied  by  changing  the  number  of  resolution  cells  averaged.  Clutter 
patch  resolution  is  determined  by  simply  multiplying  the  pixel  size  by 
the  averaging  number.  For  example,  10  x 10  averaging  will  result  in 
an  effective  resolution  of  the  clutter  patch  of  15  m by  15  m.  The 
lxl  averaging  is  an  exception,  however,  because  the  test  resolution 
is  restricted  by  the  image  resolution  of  3 meter  x 4.5  meter.  There- 
fore, the  clutter  patch  size  for  lxl  averaging  is  3 m by  4.5  m, 
not  1.5  m x 1.5  m. 

The  output  of  the  image  dissector  facility  is  9-track  800  BPI 
CCT.  Each  tape  is  given  a unique  number  that  is  recorded  in  the  ERIM 
image  processing  laboratory.  This  number  is  called  the  IPL  number. 

Every  scene  digitized  i^  recorded  in  a separate  file  on  a tape.  The 
length  of  the  file  is  proportional  to  the  size  of  the  digitized  scene. 
Each  file  is  terminated  by  a single  tape  mark.  The  first  record  of  each 
file  consists  of  120  byte  ASCII  characters.  Tlese  characters  are  used 
to  describe  the  contents  of  the  file  and  are  keyed  in  by  the  operator. 

The  data  files  are  composed  of  n records  where  n is  the  number  of  range 
scan  lines  to  be  digitized.  The  length  of  each  record  is  determined 
by  the  number  of  image  elements  of  pixels  contained  in  a single  range 
scan  line  and  each  pixel  is  represented  by  an  8-bit  word.  The  maximum 
record  length  is  restricted  to  3000  pixels  by  the  image  dissector  tube. 
The  maximum  scan  length  of  the  dissector  is  30  mm.  A schematic  diagram 
of  the  CCT  tape  format  is  shown  in  Figure  11 

2.2  DATA  PROCESSING 

2.2  1 CALIBRATION  OF  RADAR  DATA 

In  order  to  predict  performance  of  an  imaging  radar  system,  the 
designer  must  establish  values  of  numerous  radar  parameters  and  calcu- 
late signa]  levels  within  the  system.  Many  of  these,  such  as  transmitter 
power,  antenna  gain,  signal  compression  ratios,  and  noise  levels,  can 
be  determined  within  relatively  small  tolerance.  However,  target 
reflectivity  at  the  radar's  wavelength  is  a factor  that  exhibits  large 
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variation.  Airborne  radars  used  to  image  features  of  the  earth's 
surface  encounter  a variety  of  reflectivities.  Many  objects  reflect 
specularly  at  radar  wavelengths.  Some  of  these  reflect  away  from  the 
radar  and  produce  weak  echos  in  the  radar,  and  some  reflect  toward  the 
radar  and  produce  very  strong  echo  signals  in  the  radar.  There  are 
reflecting  surfaces  that  have  high  conductivity  and  reflect  nearly  all 
of  the  incident  radar  illumination  and  there  are  surfaces  that  absorb 
radiation. 

In  addition  to  the  large  dynamic  range  of  radar  echo  signals, 
there  is  uncertainty  due  to  limited  measured  values.  The  radar  cross 
section  (RCS)  of  simple  objects  can  be  calculated  from  the  size  and 
geometry  of  the  object,  the  radar  wavelength,  and  illumination  angle 
[2].  Of  course,  when  the  reflecting  object  becomes  more  complex, 
the  mathematical  model  increases  in  complexity.  Expressions  for  RCS 
require  knowledge  of  surface  roughness,  incident  angle,  complex 
dielectric  constant,  wavelength,  and  polarization.  It  has  been 
found  useful  to  consult  experimentally  derived  listings  for 
various  target  types  [3,4]. 

Although  there  is  no  built-in  calibration  in  this  radar  system, 
calibration  of  the  data  was  attempted  by  including  objects  of  known 
RCS  within  the  imagery.  All  of  the  radar  imagery  obtained  in  Labrador 
utilized  only  two  altitudes  with  their  corresponding  two  antenna  depres- 
sion angles.  An  array  of  corner  reflectors  was  installed  at  Goosebay 
Airport  and  imagery  of  the  array  was  obtained  at  each  of  the  two 
depression  angles.  A radar  image  of  the  Goosebay  corner  reflector 
array  is  shown  in  Figure  12.  This  array  was  constructed  to  provide 
calibration  points  in  the  imagery.  The  size  and  orientation  of  these 
five  corners  is  known  and  their  expected  reflecting  areas  can  be  cal- 
culated. The  reflectors  were  all  the  same  size,  but  the  illumination 
angles  differed  enough  to  provide  considerable  variation  of  effective 
radar  cross  section.  The  values  are  shown  in  Table  1. 

The  absolute  value  of  the  numbers  comprising  a digitized  image  is 
a function  of  digitizing  parameters  such  as  laser  power,  amplifier 

2.  G.  T.  Ruck,  "Radar  Cross  Section  Handbook,"  Plenum  Press,  1970. 

3.  R.  L.  Cosgriff,  W.  H.  Peake,  R.  C.  Taylor,  "Terrain  Scattering 
Properties  for  Sensor  System  Design,"  Engineering  Experiment 
Station  Bulletin  IP  Ohio  State  University,  1960. 

4.  Maurice  W.  Long,  "Radar  Reflectivity  of  Land  and  Sea,"  D.  C. 

Heath  and  Company. 


19 


*0  • 


igure  12.  Radar  Image  of  Corner  Reflectors  at  Goosebav  Airport. 


Table  1.  L-Band  Radar  Cross  Section  of  Corner 
Reflector  at  Goosebay 


Corner 

Depression  Angle 

Radar 

Cross 

Section 

A 

44.06° 

1845 

- 6.4 

dB  = 

423  sq. 

ft. 

B 

41.88° 

1845 

- 6.8 

dB  = 

385  sq. 

ft. 

C 

30.12° 

1845 

- 1.8 

dB  = 

1219  sq. 

ft. 

D 

27.46° 

1845 

- 6.2 

dB  = 

443  sq. 

ft 

E 

26.36° 

1845 

- 6.2 

dB  = 

443  sq. 

f t . 

gain,  phci.o  multiplier  tube  calibration,  and  integration  time.  Most 
of  the  digitizing  parameters  were  held  constant  for  all  of  the  digitized 
data.  However,  the  laser  power  and  the  integration  time  were  changed 
to  accomodate  high  intensity  images  in  the  city  and  in  the  corner  reflec- 
tor array  at  Goosebay  Airport.  The  recorded  numbers  were  normalized 
to  a laser  power  of  500  milliwatts  and  an  integration  time  of  900  micro- 
seconds . 

It  is  desirable  to  calibrate  these  recorded  digital  numbers  in 
terms  of  radar  cross-section  or  target  reflectivity.  Images  of  reflec- 
tors of  equal  radar  cross  section  should  vary  across  the  image  swath 
due  to  a number  of  factors.  These  factors  include: 

1.  antenna  gain  variation  with  depression  angle, 

2.  echo  signal  attenuation  with  range, 

3.  image  compression  ratio  variation  with  range,  and 

4.  diffraction  efficiency  variation  across  the  signal  film. 

The  magnitudes  of  these  factors  h ive  been  measured  and  calibration 
curves  can  be  drawn.  Table  2 lists  these  factors  at  several  points 
across  the  image  swath.  These  values  were  calculated  for  the  altitude, 
antenna  depression  angle,  and  range  used  by  the  ERIM  radar  on  February  23, 
1977,  in  Pass  5. 

A relative  calibration  curve  using  the  values  of  Table  2 is  shown 
in  Figure  13.  This  curve  indicates  that  there  is  about  14  dB  of  image 
intensity  variation  across  the  parallel  polarized  L-band  imagery  due  to 
knovm  radar  characteristics.  The  downward  slope  at  smaller  depres- 
sion angles  is  due  to  the  signal  reduction  with  range  which  has  not  been 
compensated  for  by  antenna  gain  variation.  Signal  reflection  from 
portions  of  the  aircraft  (multipath  radiation)  may  also  contribute  to 
this  signal  variation.  Diffraction  efficiency  values  were  obtained  by 
interpolation  from  measured  values. 

This  relative  calibration  data  in  conjunction  with  the  reflecting 
area  values  (Table  1)  were  used  to  compute  the  intensity  values  of  Table  3. 
The  measured  values  of  image  intensity  obtained  from  the  digitized  data 
were  also  added  to  Table  3,  and  the  differences  between  the  calculated 
and  the  measured  relative  values  shown. 
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Table  3.  Computed  and  Measured  Intensity  Values 


Corner 

Comp  Int. 

Meas.  Int. 

Meas .-Com 

A 

6.04  dB 

0 dB 

-6.04  dB 

B 

5.73  dB 

+7.33  dB 

+1.61  dB 

C 

8.84  dB 

+9.14  dB 

+0.30  dB 

D 

.86  dB 

+6.07  dB 

+5.21  dB 

E 

.26  dB 

+2.23  dB 

+1.97  dB 

If  the  reflectivities  of  the  corner  reflectors  were  as  calculated 
and  the  relative  calibration  curve  was  correct,  the  differences  between 
the  calculated  and  measured  values  should  have  been  zero.  There  was  a 
difference,  however,  and  this  difference  was  assumed  to  be  due  to  re- 
duced reflecting  area  of  individual  corner  reflectors.  The  effective 
area  of  a corner  reflector  can  only  be  decreased  from  the  theoretical 
value  due  to  mechanical  diff iciencies.  With  this  practical  considera- 
tion in  mind,  one  can  use  corner  D,  which  has  the  greatest  positive 
difference,  as  a reference  by  assuming  its  area  is  correct,  and  then  scale 
the  other  corner  reflectors  with  respect  to  D as  shown  in  Table  4. 

This  scaling  appears  reasonable  although  the  11.25  dB  reduction 
for  A is  rather  large.  This  corner  reflector  was  probably  damaged  or 
misaligned.  If  corner  D is  truly  exhibiting  its  calculated  radar  cross 
section  of  443  sq.  ft.,  an  absolute  calibration  curve  could  be  derived 
from  the  relative  calibration  curve  by  labeling  it  443  sq.  ft.  The 
computer  print-out  numbers  can  be  associated  with  this  curve  by  placing 
the  corner  D computer  print-out  number,  179,  on  the  corner  D spot  at 
27.46°  depression  angle  as  indicated  in  Figure  14.  A scale  of  numbers 
proportional  to  the  square  root  of  intensity  has  been  inserted  at  the 
right  of  Figure  14.  These  numbers  are  also  stored  on  the  digital  tape. 
Curves  for  radar  cross-sections  of  1,000  sq.  ft.,  100  sq.  ft.,  and  10 
sq.  ft.,  have  been  added  to  Figure  14. 

These  curves  can  be  used  to  determine  the  radar  cross  section  of 
any  cell  in  the  radar  image.  Interpolation  is  facilitated  by  the  logar- 
ithmic scale. 

The  information  of  Figure  14  has  been  rearranged  in  Figure  15 
to  make  the  computer  intensity  numbers  easy  to  read.  Interpolation 
in  this  figure  requires  use  of  a square  root  relationship  between 
numbers  and  image  intensity. 

For  working  with  radar  imagery  of  extended  target  surfaces 
rather  than  point  reflectors,  it  has  been  found  helpful  to  consider 
target  reflectivity  rather  than  target  cross  section.  Reflectivity 
is  intended  to  represent  a target  characteristic  without  dependence 
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Table  4 


Apparent  Reduction  in  RCS  of 
Corner  Reflectors 


Corner 

Reduction 

A 

-11.25  dB 

B 

- 3.60  dB 

C 

- 4.91  dB 

D 

0 dB 

E 

- 3.24  dB 
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upon  radar  parameters.  Reflectivity  is  defined  as  the  ratio  of  radar 
cross  section  to  target  intercept  area  [3]: 

Y = f a) 

i 

Specular  reflectors  are  observed  to  have  reflectivities  much 
greater  than  1.0  while  absorbing  surfaces  may  have  reflectivities  much 
less  than  1.0. 

Extensive  horizontal  surfaces  within  the  radar's  imaged  swath  pro- 
duce images  that  extend  over  many  resolution  cells.  In  studying  imagery 
of  cropland,  snow  covered  areas,  or  ice  packs,  reflectivity  of  each 
resolution  cell  is  of  interest.  The  reflectivity  of  horizontal 
surfaces  can  be  expressed  as: 


Y 


prpcr“"  6 


(2) 


where  P^  is  range  resolution 

PCR  is  cross  range  re  ution,  and 
9 is  depression  a Ae. 

The  depression  angle  used  in  this  expression  is  the  angle  between 
the  radar  line-of-sight  and  the  local  terrain  surface.  Terrain  that 
slopes  toward  the  radar  has  an  effective  depression  angle  larger  than 
terrain  which  is  level  or  horizontal.  The  radar  cross  section  of  a resolved 
patch  of  level  terrain  grows  smaller  as  the  depression  angle  grows  smaller  be- 
cause the  intercepted  area  of  the  incident  radiation  is  smaller.  Many  of  these 
surfaces  have  been  found,  however,  to  have  reflectivity,  y,  nearly  independent 
of  depression  angle.  In  these  cases  a single  y value  characterizes  the 
surface  of  nearly  any  depression  angle.  This  horizontal  surface  rela- 
tionship does  not  hold  at  steep  depression  angles  where  the  radar  resol- 
ution cell  is  larger  than  the  surface  intercepted  by  the  radar 
beam. 


30 


i 

! 


The  cor. -.tant  o curves  of  Figures  14  and  15  could  be  changed  to 
constant  y curves  by  the  relationship: 

a = Y PRPCR  tan  6 ('3') 

In  the  digitized  Imagery  of  Goosebay  Air  Force  Base  is  3 m and 
is  4.5  m.  The  value  of  y and  its  accompanying  computer  number  N has 
been  calculated  and  entered  in  Table  5.  Computer  numbers  associated 
with  adjacent  constant  y values  have  been  calculated  and  included  in 
Table  5. 

Figure  16  is  a plot  of  computer  numbers  for  constant  horizontal 
surface  reflectivity  over  a range  of  depression  angles.  This  figure 
can  be  used  to  find  a reflectivity  value  for  any  resolution  cell  of  tna 
radar  image  from  the  digital  numbers. 

Another  relationship  used  with  imagery  of  horizontal  surfaces  is 
reflectivity  per  resolution  cell  projected  onto  the  horizontal  plane. 
This  is: 


o 

o 


cos  0 


(4) 


Values  of  o are  a function  of  local  depression  angle  and  are 
o 

meaningless  without  identification  of  the  depression  angle  to  which 
they  apply. 


t 


Numbers  associated  with  constant  values  were  calculated  and 

entered  into  Table  6.  These  values  were  plotted  in  Figure  17  which 

shows  intensity  numbers  over  a range  of  depression  angles  for  three 

constant  a values, 
o 

The  information  of  Figure  17  was  rearranged  to  show  the  square 
root  relationship  of  intensity  numbers  with  projected  reflectivity 
values  at  a selection  of  depression  angles.  This  is  shown  in  Figure  6 
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Figure  17.  Plot  of 
Angle  f 


Figure  IB  represents  the  radar  cross  section  per  unit  ground  area, 

a , as  a function  of  the  computer  numbers.  Curves  are  drawn  for  the 
o 

three  angles  used  in  the  clutter  distribution  statistics.  Because  the 
bright  corner  reflectors  would  saturate  the  image  dissector  the  laser 
power  had  to  be  attenuated  from  500  m watts  to  50  m watts.  This  nor- 
malization is  reflected  in  Table  7 where  normalizing  factors  for  the 
digital  tapes  are  given. 

To  convert  the  computer  numbers  to  o values  you  find  the  number 

on  the  vertical  axis  of  Figure  18  and  draw  a straight  line  across  until 

you  intersect  the  curve  of  the  proper  depression  angle  and  read  the 

value  for  a off  of  the  X-axis.  This  curve  is  normalized  to  500  milli- 
o 

watts  laser  power  and  an  integration  time  of  900  microseconds. 

2.2.2  CLUTTER  DISTRIBUTION  SOFTWARE  DOCUMENTATION 

This  section  describes  the  interactive  Fortran  IV  program  which 
generates  relative  frequency  and  cumulative  distributions  from  digitized 
radar  data.  This  program  also  tests  these  frequency  distributions 
against  three  theoretical  distributions,  namely,  the  normal,  log-normal, 
and  gamma  functions  for  goodness  of  fit.  Appendix  II  lists  the  main 
program  and  associated  subroutines. 

The  basic  block  of  data  used  by  the  distribution  program  contains 
an  ASC  II  header,  which  identifies  the  data,  followed  by  the  actual 
binary  data.  Each  data  point  is  a number  proportional  to  the  square 
root  of  a radar  intensity  return  the  size  of  one  pixel.  Each  data 
block  is  a rectangular  patch  of  points  containing  typically  256  pixels 
in  range  and  4660  pixels  in  azimuth.  These  blocks,  therefore,  usually 
contain  well  over  a million  points.  In  order  to  economize  on  data 
storage,  each  point  is  first  divided  by  16  and  then  written  on  computer 
compatible  tape  (CCT)  in  8 bit,  unsigned  binary. 

The  distribution  program  is  basically  an  averaging  and  distribution 
plotting  tool  with  provisions  for  chi-squared  testing.  Figure  19  shows 
the  program  flow  which  may  be  divided  into  four  major  parts: 
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Table  7.  Normalization  Factors  For  All 
IFL  Tapes 
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SECTION  6 


INPUT  USER  RETURN  CODE 


Figure  19.  Computer  Program  Flow  Chart 
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1.  The  proper  data  is  selected. 

2.  Any  required  averaging  or  preprocessing  is  performed. 

3.  The  frequency  distribution  is  calculated  and  tested  against 
three  theoretical  curves  for  goodness  of  fit. 

A.  The  results  may  be  plotted  or  saved  for  external  post- 
processing . 

The  first  major  portion  of  the  program  includes  sections  1,  2,  and 
3 of  the  program  listing  (Appendix  II) . Initially,  the  user  is  prompted 
to  enter  the  location  of  the  data.  The  program  will  echo  the  header 
associated  with  the  data  and  then  ask  the  user  to  verify  that  the  pro- 
per data  has  been  selected.  The  data,  however,  is  not  read  at  this  time. 
The  user  must  first  indicate  whether  the  raw  data  is  to  be  squared  or 
converted  to  decxbels.  Also,  a subset  .of  the  entire  data  block  must  be 
defined  before  actual  processing  can  begin  since  a data  block  may  con- 
tain more  than  a million  points,  but  the  program  can  only  handle  a 
quarter  million  at  a time.  The  subblock  of  data  is  specified  by  entering 
azimuth  and  range  swaths.  It  is  then  that  the  program  will  actually 
read  the  specified  data  into  memory.  Although  the  data  is  physically 
stored  in  a one-dimensional  array,  it  is  organized  in  two  dimensions, 
range  and  azimuth. 

The  next  major  portion  of  the  program  is  identified  as  Section  A in 
the  program  listing.  In  this  section  the  user  will  be  prompted  to  enter 
averaging  intervals.  If  none  are  entered,  no  averaging  will  be  per- 
formed. If  intervals  are  supplied,  the  total  number  of  points  in  each 
averaging  cycle  will  be  the  product  of  the  range  and  azimuth  intervals. 
This  has  the  effect  of  decreasing  resolution  in  proportion  to  the  inter- 
val size.  For  example,  if  the  resolution  of  an  area  10  pixels  in  azi- 
muth and  10  pixels  in  range  were  to  be  halved,  square  data  groups  of 
A pixels  each  would  be  averaged  together.  This  averaging  would  reduce 
the  original  10  by  10  matrix  to  a 5 x 5 matrix.  Note  that  the  azimuth 
and  range  intervals  are  independent  of  each  other. 

After  the  data  has  been  averaged  (averaging  one  pixel  square 
sections  of  the  data  is  equivalent  to  no  averaging),  it  may  be  squared 
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or  converted  to  decibels  if  requested  by  the  user.  The  data  will  then 
be  reorganized  into  a one-dimensional  array  for  the  remainder  of  the 
program. 

The  third  major  portion  of  the  program,  identified  as  the  first 
half  of  Section  5,  will  actually  compute  the  frequency  distribution  of 
the  data  and  compare  it  with  theoretical  distributions.  Also  computed 
in  this  section  are  the  mean  and  standard  deviation  of  the  data.  These 
computations  are  performed  in  double  precision  to  ensure  no  loss  of 
accuracy  in  tLa  standard  deviation  calculation.  After  these  statistics 
have  been  displayed,  the  user  is  prompted  to  enter  the  minimum  and 
maximum  values  of  the  data  to  be  included  in  the  distribution  along  with 
the  number  of  intervals  to  be  used.  The  user  may  also  optionally 
specify  scales  for  plots  and  whether  points  lying  outside  the  min-max 
of  the  distribution  should  be  considered  equivalent  to  the  min-max 
values.  After  all  the  user  input  is  accepted,  the  program  will  calcu- 
late the  actual  frequency  distribution. 

Once  the  raw  data  distribution  has  been  calculated,  theoretical 
distributions  are  computed  given  the  mean  and  standard  deviation  of  the 
raw  data.  If  no  decibel  conversations  were  performed  on  the  raw  data, 
three  theoretical  curves  will  be  computed.  These  are  the  normal,  log- 
normal and  gamma  distributions.  If  decibel  conversions  were  performed, 
only  the  normal  distribution  will  be  computed.  These  distributions  are 
computed  by  integrating  the  area  under  the  respective  theoretical  curves. 
Due  the  the  integration  technique,  sections  of  curves  with  very  steep 
slopes  may  produce  small  negative  areas.  This  can  superimpose  small 
oscillations  on  the  cumulative  plots.  After  all  the  theoretical  distri- 
butions have  been  computed,  they  are  compared  for  goodness  of  fit  with 
the  raw  data  distributions  using  the  chi-square  technique.  The  max 
chi-square  value,  which  is  displayed  with  the  chi-square  results,  is  only 
an  approximate  figure  if  the  degrees  of  freedom  are  over  30.  All  the 
distributions  may  be  optionally  normalized  at  this  point. 
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The  last  major  portion  of  the  program,  sections  5 and  6,  will  dis- 
play the  results  of  the  previous  computations.  The  user  may  request 
a histogram  of  the  frequency  and/or  cumulative  distribution  of  all  the 
data  plotted  simultaneously.  Optionally,  a tabular  listing  of  the 
results  may  be  requested.  Finally,  the  raw  data  for  this  run  may  be 
saved  in  a file  to  be  processed  by  other  programs,  such  as,  an 
autocorrelation  program. 

The  user  is  then  given  the  opportunity  to  go  back  to  any  section 
of  the  program  for  continued  processing  after  the  cost  of  the  run  just 
completed  has  been  displayed. 

The  program  documented  in  this  section  is  a powerful,  cost-effective 
tool  for  statistical  analysis  of  radar  imagery.  It  will  process  up  to 
a quarter  million  radar  return  samples  interactively  with  the  user. 

The  histograms  produced  by  the  program  are  both  flexible  and  easy  to 
read.  Finally,  and  perhaps  most  important,  the  user  of  the  program  is 
able  to  return  again  and  again  to  any  section  of  the  program  to  modify 
or  repeat  any  portion  of  his  analysis. 


3 

STATISTICS  FOR  TERRAIN  CLUTTER 

During  the  flight  of  2/24/77  over  the  Labrador  Sea  area,  the  ERIM 
dual  frequency,  dual  polarization  synthetic  aperture  radar  measured 
clutter  return  from  several  different  kinds  of  terrain  and/or  ice 
surfaces.  This  section  makes  an  abbreviated  statistical  examination, 
at  the  L-band  frequency  for  HH  polarization  only,  of  data  recorded 
for  refrozen  lead  and  pack  ice,  black  ice,  mountainous  terrain 
(both  with  and  without  lakes),  and  tundra.  Data  recorded  for 
several  fields  and  a city  on  a separate  flight  are  also  included  for 
comparison  purposes. 

This  section  defines  the  statistics  used  to  examine  these  data 
(section  3.1);  it  also  presents  and  discusses  those  distributions 
which  are  leading  candidates  to  characterize,  on  the  basis  of  the  chi- 
square  *-.est,  the  data  histograms  (section  3.2). 

Section  3.3  contains  histograms  of  data  for  the  various  clutter 
types:  pack  ice  and  refrozen  lead,  black  ice,  mountainous  terrain, 

tundra,  fields,  and  a city  (Toronto).  In  addition,  the  effect  on  the 
data  histograms  of  averaging  over  various  sized  boxes  of  data  (or 
equivalently,  of  increasing  the  resolution  cell  size)  is  examined. 

And  finally,  in  section  3.4,  the  correlation  results  are  presented. 

3.1  STANDARD  STATISTICS  AND  CORRELATION 

The  mean 


n 


i=l 


and  stana;-:rd  deviation 
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are  computed  from  the  digitzed  data  for  each  sample  n of  data  to  be 
examined.  In  addition,  histograms  for  these  data  are  also  determined 
in  the  form  of  frequency  of  occurrence  of  n^  vs.  n^,  or  frequency  of 
occurrence  of  the  n.  in  a prescribed  interval  vs.  that  interval.  No 
higher  order  statistics  of  the  data  are  computed,  since  all  distributions 
chosen  as  candidates  to  characterize  the  histograms  are  completely 
specified  by  mean  and  standard  deviation  alone. 

To  study  the  effects  of  resolution  size  on  the  nature  of  the 

histograms,  the  data  are  averaged  for  a clutter  patch  size  determined 

\ 

by  certain  azimuth  and  range  intervals.  The  results  are  then  used 
to  form  histograms  as  indicated  above. 

Finally,  the  correlation  coefficient  p(k)  has  also  been  deter- 
mined for  various  clutter  types.  The  correlation  coefficient  is  defined 
as 


pOO 


n 


(7) 


where  n,  n.,  y,  and  0 are  as  defined  above,  and  i + k = i (modulo  n) . 
This  latter  device  has  two  effects:  (1)  it  lines  up  columns  of  range 

numbers  in  a single  line;  and  (2)  it  folds  the  last  k of  the  range  num- 
bers back  on  the  first  k of  the  range  numbers.  The  results  could  be 
slightly  misleading  when  the  terrain  corresponding  to  the  last  range 
column  is  vastly  different  from  the  first.  This  consideration  does 
not  arise  here,  however,  since  n is  always  taken  to  be  very  much  larger 
than  k;  thus,  any  slight  tendency  toward  misleading  information  is 
vastly  outweighed  by  the  programming  simplicity  required  for  the  above 
definition  of  p. 
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3.2  FITTING  DATA  WITH  STANDARD  FREQUENCY  DISTRIBUTIONS 


i 


This  section  examines  the  question  of  characterizing  the  various 
types  of  clutter  data  by  means  of  well-known  distributions.  Three 
distributions  have  been  examined  in  some  detail  in  connection  with  each 
type  of  clutter  data.  They  are  the  normal,  log  normal  [5],  and  gamma; 
each  is  defined  and  discussed  in  further  deta.il  in  section  3.2.1 
below.  These  distributions  were  chosen  for  examination  because  each 
has  appeared  in  connection  with  some  type  of  ' "rain  or  sea  clutter 
characterization  (the  Rayleigh  distribution  is  a special  case  of  the 
gamma  distribution).  There  are,  of  course,  other  distributions  which 
have  appeared  in  connection  with  clutter  studies  but  for  which  time 
limitation  has  permitted  no  study  in  relation  to  the  present  data. 

They  are  enumerated  in  section  3.2.2.  - 

Any  of  the  distributions  to  be  mentioned  in  this  discussion  can 
be  modified  somewhat  by  introducing  a "component"  with  a standard 
deviation  which  differs  from  the  standard  deviation  of  Jie  original 
distribution.  The  result  is  called  a "contamirafed  distribution." 
Although  the  term  appears  in  the  literature  [■']  with  reference  only  to 
normal  distributions,  there  is  no  reason  why  the  concept  cannot  be 
extended  to  other  distributions,  and,  in  particul.  r,  to  the  log-normal 
and  gamma  distributions.  This  extension  is  made  in  secuon  3.2.3 
for  the  gamma  distribution  and  the  resulting  improvement  of  fit  over 
an  ordinary  gamma  distribution  is  illustrated  and  discussed  for  moun- 
tainous terrain  data. 

5.  Only  about  100  years  old.  See  J.  Aitchison  and  J.  A.  C.  Brown,  The 
Log-Normal  Distribution,  Cambridge  Univ.  Press,  1957. 

6.  The  author  is  uncertain  of  the  origins  of  this  terminology  but  has 
found  it  in  connection  with  "contaminated  normal  clutter’  in  G . V. 
Trunk  and  S.  F.  George,  "Detection  of  Targets  in  Non-Gaussian  Sea 
Clutter,"  IEEE  Trans.  AES,  Vol . AES-6,  pp.  620-628,  Sept.  1970. 
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The  term  "improvement " is  used  in  connection  with  the  chi  squared 
test.  Time  has  not  permitted  the  examination  of  any  other  test,  altuough 
there  are  several  which  might  be  appropriate.  A discussion  of  the 
question  of  testing  the  fit  of  known  distributions  to  the  clutter 
data  is  given  in  section  3.2.4. 


3.2.1  NORMAL,  LOG-NORMAL,  AND  GAMMA  DISTRIBUTIONS 
In  the  present  work  the  normal  and  log-normal  distributions  are 
used  in  essentially  equivalent  senses.  The  data  are  digitized  so  that 
a numerical  value  describing  a given  pixel  of  clutter  is  proportional 
to  field  strength,  or  square  root  of  power.  Thus,  if  one  wishes  to 
determine  whether  returned  power  is  distributed  log-normally,  there 
are  two  choices  of  procedure.  On  the  one  hand,  the  digital  data  nay 
be  modified  by  applying  20  log*  to  each  number.  The  modified  data 
are  then  tested  against  a normal  distribution  whose  mean  and 

standard  deviation  a2Q  are  found  from  the  20  log  data.  Figure  20 
shows  the  normal  curve 


N(x)  = 


(x  y20  lg)  /2°20  lg 


(8) 


/2tT 


20  lg 


superposed  on  a histogram  of  pack  ice  data;  A is  the  area  under  the 
histogram.  Figure  21  is  an  alternative  presentation  of  the  same  data 
on  probability  paper.  The  staircase  represents  cumulative  frequency; 
the  straight  line  is  a cumulative  distribution  corresponding  to  Eq.  (8). 

On  the  other  hand,  a log-normal  distribution  may  be  tested 
against  the  unmodified  data  by  taking  20  log  x to  be  the  independent 
variable  of  the  log-normal  distribution.  The  resulting  log-normal 

* The  choice  of  20  log,  both  here  and  with  the  log-normal  distribution 
below,  has  been  made  so  as  to  deal  with  data  expressed  in  dB  relative 
to  some  power  level. 
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Figure  20. 


Histogram  and  Normal  Density  Fit  for  Pack  Ice 
( IPL  407,  File  1.  Pass  2,  L-Band) 


Figure  21.  Cumulative  Distribution  for  Pack  Ice 
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is  shown  in  Figure  22  superposed  on  a histogram  for  the  same  pack  ii.e 
data  used  for  Figure  20;  and  Figure  23  shows  the  same  information  on 

log-probability  paper.  The  quantities  I^q  lg  ’ a20  lg  ’ and  A ^ave  t^ie 
same  meanings  as  above;  In  denotes  the  natural  logarithm  (in  contrast 
to  log  which,  of  course,  denotes  the  common  logarithm).  Thus,  Figures 
20,  21,  and  22  are  equivalent  representations  for  the  same  information; 
the  log-normal  representation,  together  with  others,  is  chosen  for 
comparison  with  the  data  distributions- to  be  examined  in  section  3.3 
below. 


In  many  studies  a Rayleigh  distribution  has  been  used  to  fit 
certain  data.  That  distribution  is  included  here  under  the  more  general 
gamma  distribution 


g (x)  = aV  xV  1 e ax,  a,v  > 0,  x > 0 

&a,v  r(v) 


(10) 


fhe  mean,  variance,  and  mode  of  this  distribution  are,  respectively, 


\)  2 v . V - 1 

^ = a’  ° = ~V  3nd  Xm  = 
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and  the  maximum  value  of  the  gamma  frequency  distribution  at  the  mode  is 


* For  this  and  subsequent  distributions  which  are  defined  only  for 

x > 0,  it  is  always  to  be  understood  that  the  distributio  , is  defined 
to  be  0 for  x < 0. 
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Figure  22.  Histogram  and  Fitting  Log-Normal  Distribution 
for  Pack.  Ice 


2) 


The  parameters  a and  V can  be  expressed  in  terms  of  two  of  the  three 

quantities  y,CT  and/or  x . The  first  two  have  been  chosen  for  the  present 

m 

purposes:  Thus, 


a 


(13) 


and  so  Eq.  (3)  can  be  rewritten  as 
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(14) 


Now  it  is  easy  to  see  that  if  V = 1,  the  gamma  density  (14)  becomes  a 
Rayleigh  density,  with  y = CT: 


8l/a,l 
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-x/ct 
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x > 0 


(15) 


The  log-normal  distribution  (8)  and  the  gamma  distribution  (14) . 
together  with  the  standard  normal  distribution  are  used  in  section  3.3 
below  to  compare  with  histograms  of  various  types  of  Labrador  Sea  data. 
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3.2.2  OTHER  DISTRIBUTIONS 

There  are  still  other  distributions  which  have  been  used  with 
success  by  some  researchers  to  describe  certain  clutter  situations. 
Although  time  has  not  permitted  their  application  in  this  work,  they 
should  certainly  not  be  excluded  from  any  further  study. 

The  main  distributions  which  should  be  mentioned  are  the  Beckmann 

[7],  the  Rice,  and  the  Hoyt  distributions.  Briefly,  these  may  be  defined 
as  follows.  Let  rectangular  variables  x,y  and  polar  variables  r ,ip 
be  related  as  usual  and  let  x and  y be  normally  distributed  with  mean 
values  a and  3 and  variances  s and  s^,  respectively.  Then  the  dis- 
tribution for  r,  which  has  been  called  the  Beckmann  distribution,  is 


B(r) 
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where  E = m i and  I is  the  modified  Bessel  function, 

m ( 2,  m f 0,  m 

If  a = 0 and  s = s = n/2  (I  (0)  = m 7 ° ) , then  Eq . (16)  reduces 
1 2 m ( u , m f U 

2 

to  the  Rayleigh  distribution  [compare  (15)  with  x = r ] 


B(r) 
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(17) 


If  a c (a  constant),  and  = S2  - n^»  then  Eq.  (16)  becomes 
the  Rice  distribution 


R(r) 


2 -(r2+c2/n) 
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(18) 


7.  P.  Beckmann,  A.  Spizzichino,  "The  Scattering  of  Electromagnetic 
Waves  from  Rougn  Surfaces,"  MacMillan  Company,  New  York,,  1963. 
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Beckmann  calls  this  the  distribution  of  a constant  vector  c (along 
the  x-axis)  + a Rayleigh  vector.  Goldstein  [8]  has  made 
successful  application  of  this  distribution  to  terrain  data. 


Finally,  if  a = 0 and  s^  f s^ , the  result  is  the  Hoyt  distribution 


H(i ) = 


/sls2 


exp 


^ + s2)r 


4S1S2 


(19) 


Note  then  that  both  the  Rice  and  Hoyt  distributions  are  special  cases 
of  che  Beckmann  distribution. 


3.2.3  CONTAMINATED  DISTRIBUTIONS 

In  their  study  of  sea  clutter,  Trunk  and  George  [6],  employ  a "con- 
taminated normal"  distribution  with  success.  There  is,  of  course,  no 
reason  why  the  same  contamination  process  cannot  be  applied  to  any 
of  the  above  distributions.  It  has,  in  fact,  been  applied  with  prom- 
ising results,  at  least  for  a relatively  small  sample  (766  points) 
of  mountainous  terrain  data,  as  the  following  discussion  will  show. 

A "contaminated  gamma"  distribution  can  be  created  from  Eq.  (14) 
by  writing 
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x >_  0,  0 < y < 1 (20) 


8.  D.  E.  Kerr,  "The  Propagation  of  Short  Radio  Waves,"  MIT  Rad.  Lab. 
Series  #13,  McGraw-Hill,  1951. 
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The  second  term  is  the  gamma  distribution  (14)  (mean  p,  st.  dev.  O) 
reduced  by  the  factor  (1  - y)  (which  lies  between  0 and  1) . To  this  is 
added  the  first  ter'  which  is  also  a gamma  distribution  with  unchanged 
mean  p but  standard  deviation  modified  by  the  factor  K;  the  result  is 
then  reduced  by  the  factor  y.  If  0 < K < 1,  the  peak  of  the  gamma 
distribution  is  squeezed  upward,  while  its  sides  and  tail  are  squeezed 
inward  and  downward.  If  K > 1,  the  reverse  occurs. 

figures  24,  25,  and  26  show  contaminated  gamma  distributions  super- 
posed on  a histogram  of  766  points  representing  mountainous  terrain. 

In  Figure  24  the  gamma  distribution  is  uncontaminated;  that  is,  K = 1 

(or  y = 0) . The  result  is  marginally  acceptable  as  a fit.  Actually, 

2 

the  chi-square  test  yields  a statistic  of  x = 41.98;  the  test  statistic, 

2 

at  the  95%  confidence  level  and  25°  of  freedom  is  x^(25)  = 37.65. 

Thus,  according  to  the  chi-square  test,  the  gamma  distribution  is  not 
an  acceptable  fit  of  the  histogram. 

Figure  25  shows  a contaminated  gamma  distribution  with  K = 0.7  and 

y = 3/4  superposed  on  the  histogram.  The  "eye-ball  test"  suggests  that 

the  fit  is  better.  But  in  actual  fact  it  is  worse  according  to  the 

2 

chi-square  test,  since  it  produces  a statistic  of  X = 56.43. 

Figure  26  shows  a cortamiiM'  cd  g/V  <11  s rib  it  ion  w\ '.:v»  >'  » i.2, 

y = 1/4  superposed  ~n  the  histogram.  The  ’eye-ball  test"  suggests  >.r» j t 
the  fit  is  worse.  kit  the  chi-square  test  indicates  not  only  an  improve- 
ment in  fit  but  also  that  the  fit  is  chi-square  acceptable,  since  it 

2 

yields  a statistic  of  X = 35.49. 

The  choice  of  the  gamma  distribution  to  illustrate  the  contamina- 
tion process  was  a purely  arbitrary  one.  If  a more  extensive  study 
indicates  that  some  other  distribution  is  a more  likely  description  of 
ice  clutter  data,  then  the  contamination  process  should  be  applied  to 
that  distribution  to  obtain  more  precise  data  fits.  It  appears,  from 
the  brief  examination  of  the  application  of  the  process  to  the  gamma 
distribution,  that  the  contamination  parameters  K and  y could  be  chosen 
so  as  to  improve  the  distribution  fit  to  certain  parts  of  the  histogram, 
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in  particular,  say  at  the  tails.  Thus,  it  should  be  possible  to  design 
a computer  program  in  which  the  contamination  parameters  K and  y 
would  be  chosen  so  as  to  optimize  some  statistical  test,  such  as  the 
chi-square  test. 

The  idea  of  the  contaminated  distribution  is  appealing  because  it 
suggests  a possible  explanation  for  the  clutter  scattering;  that  two 
(or  more)  types  of  scatterers  are  causing  two  (or  more)  essentially 
different  deviations  about  some  mean  scattering  return. 

3.2.4  STATISTICAL  TESTS 

For  the  present  work,  only  one  test  (other  than  the  somewhat 
inexact  but  nevertheless  intuitively  helpful  "eyeball  test")  has 
been  used  to  determine  goodness-of-fit  of  certain  distributions  to 
data  histograms.  It  is  the  old  standby,  the  chi-square  test.  A 
description  of  the  test  can  be  found  in  almost  any  statistics  book,  but 

I 

it  is  useful  to  define  the  test  conditions  used  here. 

\ 

Let  f^  be  the  frequency  of  occurrence  of  data  numbers  in  the  ith 
of  N bins;  And  let  F.  be  the  value  of  the  distribution  being  tested  at 

1 2 

the  center  of\  the  ith  bin.  The  chi-square  statistic  X is  then  given  by 


The  test  hypothesis  H is  that  the  data  sample  is  from  the  distribution 

F,  i.e.,  as  N + the  data  sample  distribution  will  approach  the 

distribution  described  by  the  function  F.  The  test  statistic  used  for 

the  present  work  is  Xg5(df),  where  df  denotes  the  degrees  of  treedom 

(=  N-3  in  the  present  situation  since  the  total  of  frequencies,  the 

mean,  and  the  standard  deviations  of  the  data  sample  are  known)  of  the 

2 2 

chi-square  distribution.  If  X < XQr(df),  the  data  are  to  be  regarded 

yo  2 2 

as  consistent  with  H,  with  a 5%  chance  of  error.  If  X > X95(df),  then 
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H is  to  be  regarded  as  false  (or  else  a somewhat  rare  event  has 
occurred) . 

Because  of  the  form  of  Eq.  (21),  the  chi-square  test  places 
especially  strong  emphasis  on  the  tails  of  the  distribution.  This 
fact  is  clearly  illustrated  by  the  data  and  tested  distributions  pre- 
SGnt6d  in  Figures  23,  24,  and  25  above.  The  chi— square  test  has  been 
chosen  as  a criterion  to  be  studied  in  the  present  program  because 
the  data  at  the  tails  of  the  histograms  are  particularly  important 
(potential  false  alarm  sources).  However,  there  are  other  tests  which 
may  be  equally  good  or  better  and  which  may  not  offend  the  M?yeball 
test"  as  flagrantly  as  the  chi-square  test  does  in  Figure  25.  Other 
test  candidates  include  the  Komolgoroff- Smirnov,  the  Snedecor,  and 
the  Runs  tests. 

3.3  DISTRIBUTIONS  FOR  VARIOUS  TERRAIN  CLUTTER  TYPES  AT  VARIOUS 

RESOLUTION  SIZES 

This  section  contains  a discussion  of  the  statistical  description 
of  data  for  th  ; various  clutter  types:  pack  ice;  bare  mountain;  a 

city;  three  ri  _ds;  black  ice;  and  tundra.  The  vast  amount  of  material 
resulting  from  the  statistical  analysis  is  collected  in  Appendix  I. 

For  each  of  the  clutter  types  the  following  parametric  and  statistical 
information  is  recorded  (Appendix  I): 

1.  a definition  of  azimuth  and  range  swaths  (also 
displayed  as  a rectangular  "box") , 

2.  the  resulting  numbers  of  points  considered; 

3.  the  azimuth  and  range  averaging  intervals 
(i.e.,  resolution  interval  size;  1,1  indicates 

no  averaging),  together  with  a normalizing  option 
for  the  data  histogram; 

4.  the  minimum,  maximum,  mean,  standard  deviation, 
and  histogram  for  that  data  set; 

5.  superposed  on  each  histogram,  the  normal  (:) 
log-normal  (0),  and  gamma  (x)  distributions,  each 
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with  mean  and  standard  deviation  equal  to  the  ^ata 
mean  and  standard  deviation  (each  of  these  distri- 
butions is  uniquely  determined  by  specifying  its 
mean  and  standard  deviation);  and,  finally, 

6.  The  chi-square  test  statistic,  at  the  95%  confi- 
dence level,  for  the  appropriate  degrees  of  freedom 
(also  included)  and  the  chi-square  statistics  for 
the  normal,  log-normal,  and  gamma  distributions. 

A complete  description  of  the  program  and  its  workings  is  given  as 
Appendix  II. 

All  of  the  data  sets  to  which  the  computer  program  has  been 
applied  are  listed  in  Table  8,  together  with  the  figure  number 
under  which  each  appears  in  Appendix  I.  Given  with  each  data  set 
described  in  Table  8,  are  the  critical  parameters,  the  statistics, 
and  the  calibration  information  for  that  set;  these  include  the 
depression  angle,  azimuth  and  range  swaths,  number  of  points, 
resolution  size  (i.e.,  averaging  interval  size),  mean,  standard  devi- 
ation, the  chi-square  test  statistic,  and  the  chi-square  statistics 
for  the  normal,  log-normal,  and  gamma  fitting  distributions.  In 
addition,  the  calibration  information  is  used  to  convert  mean  and  stan- 
dard deviation  into  O values,  i.e.,  RCS  per  unit  illuminated  area. 

o 

Thus,  Table  8 can  be  used  not  only  as  a key  to  the  available  data  given 
in  Appendix  I but  also  as  a comprehensive  summary  of  available  results. 

The  data  figures  are  arranged  according  to  the  patterns  shown  in 
Tables  9a  and  9b.  The  arrangement  is  built  on  a basic  data  set 
which  is  comprised  of  four  items,  each  at  three  elevation  angles  for 
the  four  terrain  types:  a total  of  48  in  the  basic  set.  The  four  items 
treated  are:  a large  sample  with  no  averaging;  the  original  sample 

with  averaging  over  squares  10  pixels  on  a side;  the  original  sample 
with  averaging  over  squares  25  pixels  on  a side;  and  the  original 
sample  with  averaging  over  squares  50  pixels  on  a side.  In  addition, 
there  is,  for  the  mountain,  ice,  city  and  field  terrains,  a supplementary 
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Table  8.  Data  Inventory  (Continued) 
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Table  8.  Data  Inventory  (Continued) 
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Table  9.  Figure  Numbers 


a.  Figure 

Numbers 

for  Basic 

Data  Set 

Elevation 

Mountain 

Ice 

City 

Field 

20° 

1-4 

13-16 

25-28 

A:  37-40 

30° 

5-8 

17-20 

29-32 

B:  41-44 

60° 

9-12 

21-24 

33-36 

C:  45-48 

Each  set  of  four  figures  is  arranged  according  to  averaging  over 
squares  of  pixel  size:  lxl,  10  x 10,  25  x 25,  and  50  x 50. 


b.  Figure  Numbers  for  Supplementary  Data  Set 


Elevation 

Mountain 

Ice 

City 

Field 

20° 

49-5/ 

A: 

97-101 

30° 

53-65 

66-88 

89-96 

B: 

102-106 

05 

O 

O 

C: 

107-111 

c.  Figure  Numbers  for  Additional  Data  Set 
Elevation  Black  Ice  Tundra 


20° 

30°  112-117  118-124 

60° 
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data  set  which  is  based  on  a variation  of  sample  size  and  further 
averaging  square  sizes.  Finally,  additional  data  have  been  processed 
for  black  ice  and  tundra  terrains;  these  are  enumerated  in  part  c of 
Table  9. 

Figures  3,  4,  and  6 show,  respectively,  the  images  of  the 
ice,  mountain,  and  city  areas  which  have  produced  the  clutter  data  for 
this  study.  In  addition,  three  fields,  chosen  for  their  homogeneity 
of  image  and  medium  brightness  of  scatter,  have  been  included  for  com- 
parison purposes.  They  are  shown  in  Figures  7,  8,  and  9. 

The  depression  angles  appropriate  to  their  data  are:  Field  Af  Field  B, 

Field  C. 

The  relative  sizes  of  mean  return  and  especially  of  deviation 
about  the  means  are  shown  in  Figure  27  - The  reader  should  take  account 
of  the  fact  that  the  numbers  for  the  city  scattering  should  be  increased 
by  a factor  of  4.74  (as  is  indicated  on  the  figure)  (note  that  the  standard 
deviation  must  also  be  increased  by  the  factor  4.74).  The  mean  of 
city  data  is  relatively  insensitive  to  elevation  angle;  moreover, 
the  standard  deviation  of  city  data  is  not  only  relatively  insensitive 
to  elevation  angle  but  also  to  averaging  square  size. 

For  the  bare  mountain  and  pack  ice  data  the  mean  and  standard 
deviations  for  20°  elevation  are  quite  similar  in  behavior.  But  for  the 
other  elevations  they  are  not.  In  fact,  the  behavior  of  the  bare 
mountain  and  pack  ice  data  appears  to  be  reversed  at  30  and  60 
elevations.  More  precisely,  for  moderate  (25  x 25)  to  large  (50  x 50) 
averaging  square  sizes  the  standard  deviation  behavior  of  homogeneous 
ice  at  30°  elevation  resembles  that  of  the  bare  mountain  at  60°.  One 
can  readily  surmise  (see  next  paragraph)  that  the  bare  mountain  topo- 
graphy may  be  at  least  part  of  the  cause  of  this  interesting  standard 
deviation  behavior  with  averaging  square  size. 

The  mean  and  standard  deviation  of  the  bare  mountain  data  at  30° 
degression  angle  very  much  resemble  those  respective  statistics  for 
black  ice  at  the  same  depression  angle.  In  fact,  their  entire 
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Figure  27.  Standard  Deviation  Vs.  Resolution 


Figure  27.  Standard  Deviation  Vs.  Resolution  (Continued) 


distributions  are  quite  similar,  as  Figure  28  shows.  Now  even  though 

the  mountain  terrain  is  quite  bare  and  in  many  places  covered  by  ice 

(from  ground  truth  recorded  by  the  experimenters),  une  would  expect 

somewhat  higher  values  of  O than  those  observed  in  Figures  5 through  8 

(and  Figures  53  through  65  in  Appendix  I as  well).  Thus,  a re-examination 

of  the  mountain  terrain  was  made  (see  Figure  3) . The  area  chosen  for 

analysis,  although  quite  homogeneous  and  moderately  bright  in  the 

imagery,  appears  to  be  veiy  possibly  situated  on  a reverse  snope.  New 

sample  areas  were  chosen  for  comparison  and  these  produced  much  higher 

(10  or  more.  dB  in  0 ) mean  values.  The  resulting  wide  variation 

o 

in  average  values  from  bare  (and  somewhat  icy)  mountainous  terrain  is 
an  important  fact  to  be  noted. 

In  general,  with  the  exception  of- the  20°  elevation  cases,  the 
standard  deviation  about  the  mean  return  declines  steeply  as  averaging 
square  size  increases,  and  this  trend  is,  of  course,  to  be  expected 
since  individual  scatterers  have  less  opportunity  to  make  their 
individual  returns  stand  out. 

As  was  remarked  at  the  outset  of  this  section,  the  distributions, 
both  density  and  cumulative,  for  the  data  sets  enumerated  above  in 
Tables  8 and  9 have  superposed  on  them  a normal,  log-normal,  and  a 
gamma  distribution.  These  three  analytic  distributions  are  chosen  to 
have  mean  and  standard  deviation  equal  to  the  mean  and  standard  devi- 
ation for  the  data  histogram,  whether  or  not  any  of  the  three  distri- 
butions is  an  acceptable  fit  by  the  chi-square  test.  In  fact,  some  of 
the  chi-square  statistics  for  the  normal,  log-normal,  and  gamma  distri- 
butions are  ridiculously  high  in  comparison  to  the  test  chi-square 
statistic  for  the  data  in  question.  But  they  are  included,  whether 
they  are  high  or  not,  to  indicate  which  analytic  distribution  best 
fits  the  data,  however  badly  or  well. 

To  indicate  at  a glance  which  distributions  best  fit  the  various 
data  sets,  a tabulation  is  given  in  Table  10.  In  more  than  half  of 
the  cases,  the  log-normal  distribution  provides  the  best  fit,  with  the 
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gamma  and  normal  distributions  appearing  more  frequently  (the  gamma 
twice  as  often  as  the  normal)  as  resolution  (r,r  averaging)  size  is 
increased . 

it  is  of  considerable  interest  to  compare  not  only  means  and  stan- 
dard deviations  of  the  data,  as  has  been  done  above,  but  also  to  compare 
the  histograms  themselves  for  the  various  terrain  types.  Such  a com- 
parison is  made  in  Figure  29  for  bare  mountain,  pack  ice,  field  B, 
and  tundra  [the  data  are  actually  taken  from  Figures  54  (mountain), 

70  (pack  ice),  193  (field  B) , and  123  (tundra)  of  Appendix  I].  The 
pack  ice  distribution  has  the  longest  tail,  while  the  bare  mountain 
has  the  shortest.  Field  B and  the  tundra  distributions  are  quite 
similar  in  shape  but  the  field  B distribution  is  displaced  consider- 
ably to  the  right,  because  of  its  larger  mean,  of  the  tundra  distribution 

3.4  CORRELATION  RESULTS 

Let  ( be  a random  variable  whose  behavior  is  described  by  the 
density  function  p(£).  The  expectation  E of  a function  f(S)  is  taken 
as 


E [ f (x)  ] = J f(Q  P(D  <1S,  (22) 


where  the  integration  is  taken  over  the  defined  domain  of  S.  The 
mean  p and  standard  deviation  a used  in  this  discussion  are  well 
known  f 9 ] : 


p = E[£]  = J S p(D  d£ 

(23) 

02  = E[(S  - p)2]  = f (€  - P)2  p(x)dx 

(24) 

Let  r|  be  a second  random  variable  with  density  9(u)  . 
the  joint  central  moment 

If  pr  is 

So 

p^  = El(£  - Pj-) (n  - Pn)3, 

(25) 

9.  d.  Cramer,  "Mathematical  Methods  of  Statistics,"  Princeton 
University  Press,  Princetor,  1958. 
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Comparison  of  Histograms 


then  the  correlation  p between  the  random  variables  £ and  y is  defined 
to  be 


P = 


__S. 

a 

K n 


(26) 


where  the  subscript  denotes  the  variable  to  which  the  statistic  per-' 
tains.  If  the  expectation  y^  is  expanded,  one  obtains  the  alter- 
native correlation  formula 


P = 


E[£nl  - y^.  y^ 

or  a 

5 n 


(27) 


The  correlation  coefficient  p has  the  range  - 1 £ p £ 1. 

The  definition  (26)  affords 


y - M * - V, 

L = p ‘ 

“ r* 


(28, 


as  a best  linear  mean  square  regression  of  y by  means  of  £,  i.e., 
the  square  of  y minus  a linear  function  of  £,  viz. 


n - 


(.  ■ y) 


p O 


+ 


Oc, 


is  as  small  as  possible.  Alternatively,  a best  linear  mean  square 
regression  of  £]  by  y produces 
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. . 


. 'V  . 


(29) 


y 


l 

p 


o 


€ 


If  p is  1,  or  -1,  then  the  whole  mass  of  the  distribution  is  located 
on  a straight  line.  If  p is  0 then  (28)  reduces  to  y = ^ and  (29)  reduces 
to  x = Uj.  which  are  lines  parallel  to  the  y-  and  x-axes,  respectively, 
a.id  the  random  variables  and  n a.re  independent. 


The  above  definitions  and  interpretation  of  correlation  p have 
been  applied  to  the  clutter  data  as  follows.  Suppose  there  are  N lines 
of  data  between  two  desired  azimuths,  each  line  containing  M pixels 
corresponding  to  the  same  fixed  range  interval.  A sample  of  MN  values 
from  a random  variable  £ is  created  by  arranging  in  tandem  the  N azi- 
muth lines  (of  M pixels  each).  If  a circular  arrangement  of  the  MN 
pixels  is  made  by  closing  the  last  with  the  first,  then  the  second 
random  variable  is  generated  by  shifting  the  data  through  some  desired 
number  of  points.  The  result  is  actually  an  auto-correlation.  More 
precisely,  let  x(m),  m = 1,  2.  ...,  MN,  represent  the  sample  points, 

with  x(MN  + m)  = x(m) . Then  the  correlation  coefficient  p becomes 
(from  Eq. (27) 


p(n) 


EU(£  + n)  ] - 

2 

Or 


(30) 


The  bare  mountain  terrain  data  have  been  examined  for  two  cases 
at  the  30°  depression  angle  (both  from  IPL  411,  File  1,  Fit.  2/23/77, 
Pass  4 (L-band) . The  first  set  is  a sample  of  256  range  points  from 
azimuth  line  3301.  The  autocorrelation  coefficient 

p(n),  n = 1,  2,  . . . , 10  (31) 

is  determined  and  plotted  to  show  the  actual  values  of  p(n) . Figure 
30  contains  the  results.  For  example,  one  reads  that  p(5)  = 0.27967, 
while  p(l)  = 0.83421. 
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The  second  set  is  a sample  of  25&00  points  from  the  100  azimuth 
lines  3301  through  3400  (i.e.,  the  range  swath  is  1-256  and  the  azimuth 
swath  3301-3400).  The  results  for  p(n),  again  for  n = 1,  ....  10,  are 
shown  in  Figure  31.  This  time  one  finds  p (5)  = 0.325C'C)  and  p(l)  - 
0.85329. 
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Figure  31.  Correlation  Coefficients  p (n) , n = 1,  . ..,10  for  100 
Lines  of  Bare  Mountain  Data  (1PL  407,  File  3,  Pass  2) 
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SIMULATION  OF  RADAR  SIGNALS  FROM  AN 
R,  R CLUTTER  MAP 

The  clutter  data  statistics  were  generated  for  the  purpose  of 
generating  clutter  models  for  northern  region  terrain  types,  as  an 
input  for  the  generation  of  space  radar  clutter  in  the  form  of  mag- 
netic data  tapes.  ' A necessary  part  of  this  effort  was  to  prepare  a 
space  radar  system  model  which  simulates  the  impact  of  sensor  system 
parameters  on  clutter.  This  simulation  addressed  the  effects  of 
platform  orbit,  antenna  pattern,  waveform  bandwidth,  and  viewing  geom- 
etry. This  simulation  of  radar  signals  from  a range  (R)  and  range 
rate  (R)  clutter  map  is  described  in  this  section. 

The  procedures  developed  were  integrated  into  a computer  program 
which  was  used  to  model  the  effect  of  the  radar  transfer  function  and 
to  synthesize  time-dependent  radar  returns  from  ground  clutter  in 
order  to  create  data  tapes  which  can  simulate  signals  at  the  output 
of  the  analog-to-digital  converter.  These  tapes  contained  approxi- 
mately the  equivalent  of  150  msec  of  radar  signal  at  10  MHz  simulated 
bandwidth.  Both  in-phase  (I)  and  quadrature  phase  (Q)  channels  were 
provided . Also,  target  signature  data  was  analyzed  and  target  signals 
synthesized  in  a compatible  format  to  be  used  with  the  clutter  tapes. 

The  clutter  signal  as  a function  of  time  (t)  is  given  as 


c(t)  = co  f . MR,R)e4TrjtR/X  s(t  - 2R/c)dRdR 

A.R  (1 
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where  s(t)  is  the  waveform,  cQ  is  a complex  constant,  x is  the  wave- 
length, R is  the  range,  c is  the  speed  of  light,  R is  the  range  rate, 
j = V-]  , and  A(R,R)  is  the  complex  reflection  coefficient  of  clutter 
whose  expected  value  E|A(R,R)|^  is  the  clutter  RCS  per  unit  range  and  range  rate. 
The  latter,  the  clutter  RCS  density,  or  "clutter  map"  is  a subject 
in  itself  and  requires  careful  modeling  to  faithfully  include  effects 
of  viewing  geometry  and  terrain  type.  It  requires  antenne  patterns, 
platform  and  earth  motion,  and  viewing  angles.  Equally  important  ^nd 
least  well-defined  is  the  sigma-zero  model  which  is  also  required. 

Essentially  the  radar  cross-section  per  unit  area  for  every  point  on 
the  earth,  could  be  handled  by  maps  based  on  actual  clutter  data,  a 
stochastic  model  based  on  probability  distributions  and  spatial 
correlations  of  actual  data,  or  a suitable  combination  of  the  two. 

Realistically,  the  last  solution  is  probably  the  most  suitable  in  most 
cases. 


For  notational  simplicity  we  define  the  following  quantities 


t' 

v' 

W 


2 R/c 

2 R/X 
e-2nj 
e“27rj/N 


Hereafter,  j can  be  used  as  an  index  as  well  as  -^-T  since  its  meaning 
should  be  clear  from  context. 


For  the  moment  consider  the  waveform  to  be  a weighted  train  of 
<5  function  pulses 


N-l 

s(t)  = s0  £ wn  6(t  - nA) 
n=0 


(2) 
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where  Sg  is  a constant,  is  a set  of  complex  weights  and  a is  the  inter- 
pulse period.  The  6-function  is  defined  such  that  for  an  arbitrary  function 


co 

f g(t)  6 ( t - t')dt  = g(t') 


Substituting  (2)  into  (1)  and  integrating  yields 
c(  -)  - cnsnV;  /dv*  A(t  - nA,  v' )W-tv'w 


Now  suppose  that  c(t)  is  sampled  at  discrete  times  tg  + mA  for  N 
samples  (pulses). 

c(tg  + mA)  = CgSg^  f “dv'  A(t.g  + in  - n)A,  v' ) w",nAv ' wn 


CpSgF, j dv  A(tg  + nA,  v) W 

n -® 

= S cn(t0  + 

n 

re  cp(mA)  = CgS^dv'  A{nA,  v')WmAv  . 


(N  - 1)  < n < N - 1 


Hereafter,  the  reference  time  tQ  will  not  be  noted  unless  referred  to 
specif! cal ly. 

Consider  the  discrete  Fourier  transform  (DF7)  of  cn(mA)  given  by 
N-l 

c (v)  ~ / . \ mAv 

nK  1 Z-i  c (n»A)w 

m=0  n 

- C0S0 P"  »')E  w“('"V')  Vn 

m=0 

/«*? 

dv'  A(nA,  v')  xn  (v  - v') 


where 


N-l  . ,N 

( i\  _ V1  , .mA ( v - v ) 

x (v  - v ) = > , w w ' 

" mO  m-n 

= DFT(Vn)*  (6) 

n ~ 0,  + 1 , ...»  +^N  - 1). 

From  (5)  and  the  "additive"  property  of  the  DFT,  we  note  that 


C(v) 

= DFT(c(t) ) = 

51  cJv) 

n n 

Thus 

a procedure  for 

computing  c(t) 

is 

1. 

compute  A(R,R) 

2. 

compute  xn(^) 

- (N  - 

1 ) < n < N - 1 

by  Eqn  6 

3. 

compute  cp(v) 

- (N  - 

1)  < n < N - 1 

by  Eqn, 5 

4. 

compute  c(v) 

by  Eqn  7 

5. 

inverse  DFT  to  compute  c(t). 

This  procedure  is  undesirable,  however,  for  several  reasons  which  can 
be  remedied  by  considering  the  following  alternate  procedure. 


This  alternate  procedure  makes  use  of  the  fact  that  temporal  and 

spectral  clutter  samples  are  independent.  Unlike  the  first  method, 

this  allows  the  major  portion  of  the  work  to  be  performed  with  rear 
• 2 

quantities  a(R,  R)  and  | xn ( v ) I s where 

a(R,R)  = |A(R,R)|2  (S) 

The  bar  over  a qunatity  denotes  the  expected  yalu*  o;  a quantity, 

* Note  that  DFT  is  used  both  as  an  abbreviation  toi  i^s.-rete  Fr 

transform  and  as  a function  representing  the  if  the  • . - '.'be 

Fourier  transform  of  the  variable(s). 
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p 

In  addition  to  being  real,  \ and  c are  now  deterministic  and, 
therefore,  far  fewer  random  numbers  need  to  be  generated  and  these  only 
in  the  later  stages  of  the  process. 

Consider  the  expected  value  of  the  square  of  cn(v) 

|cn(v)]2  = lc0l2  lsol2  E[/W  A(na,  v')  A(nA,  v") 

x Xn^-V)  xj  (v  - v")|  (9) 

The  statement  that  spectral  samples  are  independent  can  be 
expressed  as 


E [A(nA,  v1)  A*(mA,  v")j  = a(nA,  v')  6(v  - v")  6fnn 

1 if  m = n 
0 otherwise 

Therefore, 

oo 

|cn(v)|2  = |c0s0|  y*dv'a(nA,  v')  X^(v  - v') 

— 00 

and  (c(v)  | 2 = ^|cn(v)|2 
n 


where 


mn 


(10) 


Ola) 

(11b) 


2 

Note  further  that  since  xn  is  computed  from  a regular,  discrete 
sample,  it  is  a periodic  function  whose  period  is  eoual  to  the  pulse 
repetition  frequency,  PRF  = f = 1/A.  Thus  it  is  possible  to  break-up 
<j(nA,  v')  into  a sum  as  follows: 
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(12) 


k = 0,  + 1 , ...+® 


where 

\ 

'k(nA,  v‘ ) = o(nA,  v'  + kfr) 

0 

1 A 
<\ 

A 

= 0 

otherwise 

Now  we  have 

f 

lCn  * IVO' 

•'a 

(nA,  v')  (v  - 

k n 

0 


(13) 


(14) 


The  right-hand-side  of  (14)  is  a periodic  convolution  and  is  most 
efficiently  done  using  the  convolution  therorem  and  fast  Fourier 
transform  (FFT).  The  convolution-theorem  is  written  as  follows: 

I. 

Zn  ” \ ^n-m  (15a) 

m=0 


where  y 
then 

is  periodic,  yp  = yn+kL> 

k = 0,  +_  1 , 

z = -r-'Y'  AmBmwi"nm 
n L m m L 

m 

= inverse  DFT(AB) 

(15b) 

where 

\ =?vC 

= DFT(x) 

(15c) 

= DFT(y) 

(15d) 

£ 
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We  apply  this  theprem  to  Eqn  (14)  with  the  approximation 


(16a) 


and  V „ mf  /L 
r 


(16a) 


The  integer  L must  be  chosen  such  that  each  "lobe"  of  xn(v)  is  sampled 
at  least  twice.  This  requires,  if  there  are  N pulses,  that  L >_  2N. 


v - 


£- 

Now  we  have  the  following  procedure  for  computing  |c(v)|  for 
m = 0,  1,  ...L-l. 


1. 


Compute  cr(R,  R) 


- 1 - n < N - 1 
0 < m < 2N  - 1 


4N 


2.  Compute  xn(v)>  n = 0,  + 1 + (N  - 1 ) 

for  v = mfr/2N,  m = 0,  ...»  2N  - 1 

p 

3.  Convolve  a with  x using  (15)  and  (16). 

4.  Sum  over  n to  obtain 


I c(v  - mfr/2N)j2  = |cj 

5.  Smooth  to  obtain  N distinct  samples  from  the  2N  available  points. 

6.  Compute  cm  = |cm|e1#1.  He*-e  |cmj  is  a Rayleigh  distributed 

random  variable  with  expected  value  Vi  cml  and  a 

random  variable  uniformly  distributed  between  0 and  2n. 
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7.  Inverse  DFT  to  obtain  the  sequence 

t 

8.  Repeat  for  every  required  value  of  tg. 

At  this  point  it  is  worthwhile  to  describe  the  physical  significance 
of  Steps  1 through  8 above. 

Step  1 is  the  production  of  a valid  clutter  map  taking  into  account 
terrain  features  and  viewing  geometry.  The  procedure  is  documented 
elsewhere  and  involves  a translation  of  area  in  earth  coordinates 
into  radar  coordinates  { R , R)  weighted  by  appropriate  factors  in 
the  radar  equation  (e.g.,  antenna  pattern,  o°  model,  and  R factor). 
For  a realistic  simulation  the  a0  model  should  reflect  a grazing 
angle  dependence  and  appropriate  terrain  mix.  In  addition, 
fluctuation  from  mean  value  and  spatial  correlations  should  be 
included  either  by  means  of  a stochastic  model  or  by  using  actual 
unsmoothed  a0  data. 

2 

Step  2 is  the  computation  of  the  cross  ambiguity  function  x 
of  the  transmitted  waveform  with  no  receive  weights  (receive 
weights  will  be  included  in  simulated  processing  of  the  simulated 
signal  produced  here).  This  is  done  by  the  FFT  algorithm. 

Step  3 is  the  convolution  of  x2  with  the  expected  clutter  in 
frequency  space.  Physically  each  point  of  the  convolution  is  the 
expected  clutter  at  that  frequency  taking  into  account  Doppler 
ambiguities  and  Doppler  sidelobes. 

Step  4 is  the  summation  of  each  of  the  ambiguous  range  strips 
into  one  equivalent  Doppler  spectrum. 

Step  5 is  simple  combining  points  which  lie  in  the  same  Doppler 
fi Iter. 

Step  6 is  the  production  of  a signal  spectrum  having  Rayleigh 
statistics  and  whose  expected  magnitude  squared  equals  that 
computed  in  Steps  1 through  5. 

Step  7 is  simply  the  computation  by  DFT  of  the  time  sequence  of 
radar  pulses  at  the  receiver. 
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Step  8 is  not  really  a separate  step  but  a reminder  that  the 
process  must  be  repeated  for  every  range  resolution  cell  to  be 
processed.  The  results  for  the  numerous  range  strips  are  not 
independent  since  they  all  come  from  a clutter  map  (Step  1)  which 
takes  into  account  spatial  correlations. 

There  is  another  step  which  could  be  called  Step  9 which  has 
not  been  mentioned  heretofore.  Since  the  pulses  have  a finite  width 
and  some  frequency  modulation,  the  single  pulse  waveform  (assuming 
each  pulse  is  identical  except  for  a complex  weight)  must  be  convolved 
with  the  time  sequence  generated  in  Steps  1 through  8.  This  is  derived 
mathematically  below.  In  essence  it  folds  together  results  of  nearby 
range  cells,  leaving  it  up  to  the  pulse  compression  process  of  the 
processing  simulation  to  untangle  them. 

This  procedure  is  valid  only  for  6-function  waveforms  of  the  type 
(3).  In  general  the  waveform  will  be  a train  of  identical  pulses  of 
finite  length.  This  is  given  by 

s(t)  = {"<■ it"  s0:t")  nt  - 1") 


where  T(t)  is  a train  of  6-function  as  in  (3),  and  s0(t)  is  the 
envelope  of  the  individual  pulses. 


Since  the  process'of  convolving  a long  sequence 
cl<t0  + l"4)  ‘min  i i ‘max’  m = 0 N ‘ 1 

1 ' 

with  a short  one  such  as  sQ{t)  (the  individual  pulse  length  is  certainly 
less  than  a and  probably  much  shorter)  is  straightforward,  it.  seems 
appropriate  to  first  compute  c^(t)  and  then  convolve  it  with  Sg(t)  by 
means  similar  to  the  FFT  convolution  algorithm  (15). 


The  above  sections  of  Volume  I and  the  appendices  in  Volume  II 
contain  Labrador  Sea  clutter  data  taken  at  L-band , HH  polarization. 

These  data  cover  pack  ice  (and  refrozen  lead),  bare  mountains,  tundra, 
and  black  ice.  For  comparison  purposes  data  are  also  included  for  a 
city  (Toronto)  and  three  homogeneous  fields.  Methods  of  measurement, 
recording,  and  processing  (including  calibration,  digitization,  and 
imagery)  of  the  data,  as  well  as  the  test  sites.,  are  described  in 
Section  2.  These  clutter  data  are,  of  course,  best  described  by  sta- 
tistical quantities.  Section  3 gives  a preliminary  examination  of 
some  of  the  statistics  by  which  the  data  can  be  characterized.  If 
the  data  are  to  be  used  to  formulate  a .mathematical  model  for  target 
detection  and/or  tracking  in  the  presence  of  clutter,  then  it  is  use- 
ful to  learn  whether  various  types  of  clutter  data  can  be  fitted  by 
analytic  distributions.  Section  3 also  addresses  this  question  and 
finds  that  the  log-normal  distribution  seems  to  be  the  number  one 
candidate  for  most  acceptable  fitting  distribution.  There  seems  to  be 
some  promise  that  an  analytic  distribution,  such  as  the  log-normal, 
can  be  contaminated  (see  section'  3 . 2 . 3)  so  as  to  produce  a distribution 
whose  fit  is  acceptable  to  a statistical  test  such  as  the  chi-square. 

The  numerical  values  of  the  data,  resulting  from  digitization  of 
recorded  signals,  are  proportional  to  the  square  root  of  intensity  (or 
received  power) . These  numerical  values  provide  a convenient  basis 
in  terms  of  which  to  do  further  processing  and  analysis.  This  approach 
has  been  used  exclusively  in  the  above  discussions.  It  is,  however, 
straightforward,  with  the  help  of  the  calibration  information,  to 
convert  these  square-root-of-intensity  numbers  into  radar  cross  section 
per  unit  illuminated  area  a . This  conversion  is  done  for  averages  of  square- 
root-of-intensity  numbers  for  pack  ice,  bare  mountain,  field,  and  city 
(Toronto)  data  in  order  that  they  may  be  compared  with  existing  data. 

Figure  32  contains  the  results,  showing  C vs.  depression  angle 
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0.  The  city  data  are  highest  valued  at  the  20°  depression  angle, 
being  aoout  -15  dB  on  the  average.  This  level  of  then  falls  off  to 
approximately  -20  dB  at  6 = 30°,  and  still  further  to  about  -24  dB  at 
0 = 60°.  It  seems  quite  reasonable  that  the  largest  values  should 
occur  at  smallest  depression  angles  since  there  the  buildings  form 
dihedrals  with  the  longest  base  sides.  However,  the  city  Gq  levels 
recorded  in  Figure  32  are  some  6 to  15  dB  below  levels  for  Phoenix 
and  a New  Jersey  residential  area  reported  by  Long  [4],  with  the  largest 
difference  occurring  at  0 = 60°  for  the  New  Jersey  residential  area. 
Rooftops  may  be  largely  responsible  for  the  higher  return  from  this 
New  Jersey  residential  area. 


The  bare  mountain  data  agree  rather  well  with  the  Arizona  mountain 
data  from  Long,  the  former  being  only  some  6 to  8 dB  lower.  Such  lower 
values  of  the  mountains  of  the  Labrador  Sea  area  could  very  well 
be  attributable  to  their  smoothness  and  to  the  fact  that  large  shadow 
areas  (on  reverse  slopes)  can  be  present  in  a given  scene. 

L-band  clutter  data  are  in  short  supply.  So  to  provide  added 
material  for  comparison,  one  sea  data  curve  from  MacDonald,  1956,  [10]  is  in- 
cluded in  Figure  32.  The  data  are  for  a sea  under  a 30  knot  wind  and 
thus  may  be  reasonably  comparable  with  the  pack  ice  data.  One  sees 
in  fact  only  a 5 or  so  dB  difference  between  them. 

The  field  data  provide  only  one  value  at  0 = 30°.  The  field  was 
homogeneous  and  was  chosen  for  its  fairly  bright  image.  Thus,  it 
should  and  does  compare  well  with  the  pack  ice  which  appeared  to  have 
about  the  same  level  of  return  in  the  imagery. 


10.  F.  C.  MacDonald,  "The  Correlation  of  Radar  Sea  Clutter  on  Vertical  and 
Horizontal  Polarizations  with  Wave  Height  and  Slope,"  1956  IRE 
Convention  Record,  Pt.  1,  pp.  9-32,  1956. 
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In  summary  then,  the  clutter  data  for  the  pack  ice,  bare  mountain, 
city,  and  field  seem  to  provide  quite  believable  O'  values  which  are  in 
reasonable  agreement  with  the  small  amount  of  L-band  clutter  data  pres- 
ently existing  in  the  literature. 

The  ERIM  clutter  data  presented  in  Figure  32  has  some  interesting 
aspects.  One  would  normally  expect  the  curves  to  decrease  with 
smaller  depression  angles.  The  city  data  exhibits  the  opposite  de- 
pendence, i.e.,  the  o’  increases  with  decreasing  depression  angle. 

This  is  probably  explained  by  the  large  flat  sided  buildings  found 
in  cities.  The  radar  return  becomes  more  specular  at  low  depression 
angles.  A similar  dependence  is  observed  on  the  pack  ice  data.  The 
pack  ice  a changes  by  approximately  3 dB  between  20  and  60  depres- 
sion angle  and  shows  a larger  a at  20°  than  it  does  at  60  . This 
would  seem  to  go  against  ones  intuition  even  though  the  difference  is 
only  3 dB.  This  small  reversal  in  the  slope  of  the  curve  could  be 
due  to  calibration  errors. 

One  of  the  fundamental  assumptions  inherent  in  the  calibration 
procedure  was  that  the  antenna  pattern  was  reasonably  accurate  and 
that  any  deviation  of  the  measured  values  of  RCS  from  the  computed 
values  for  the  corner  reflectors  was  due  to  errors  in  the  RCS  of  the 
corner  reflectors.  In  retrospect,  this  may  not  be  the  case.  The 
difference  could  have  been  due  to  antenna  gain  variations  in  the 
elevation  plane.  Again  using  hindsight,  it  probably  would  have  been 
better  to  assume  that  the  corner  reflectors  were  of  the  correct 
RCS  and  that  any  differences  between  calculations  and  measurements 
were  due  to  deviations  from  the  antenna  patterns  used  for  the  calculations. 

Section  4 contains  a description  of  a radar  signal  simulation  by 
Philip  Tomlinson  of  Decision-Science  Applications,  Inc.  This  simula- 
tion was  used  to  provide  magnetic  tape  data  for  use  in  testing  signal 
processor  techniques.  This  simulation  is  capable  of  accepting  as  input 
the  clutter  data  and  statistics  generated  on  this  program.  However, 
much  more  can  be  done  with  the  more  extensive  data  now  available. 
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CONCLUSIONS  AND  RECOMMENDATIONS 

The  data  presented  and  discussed  in  this  report  are  valuable 
additions  to  the  very  small  amount  of  existing  data  on  L-band  terrain 
clutter.  They  represent  a rather  small  percentage  of  the  total  amount 
of  data  recorded  over  the  entire  Labrador  Sea  measurement  flight. 

Thus,  there  is  available  a virtually  untapped  mine  of  data  for  various 
types  of  terrain  (and  sea)  clutter.  A sufficient  amount  of  analysis  is 
now  in-hand  to  offer  a rather  tantalizing  possibility  for  the  charac- 
terization of  these  data. 

Although  this  data  is  not  fully  amplitude  calibrated,  the  partial 
calibration  used  allows  one  to  draw  some  conclusions  with  respect  to 
this  L-band  clutter  data.  The  principle  conclusions  are: 

1.  The  clutter  statistics  generated  on  this  program  appear 
to  be  well  behaved.  This  is  a good  indication  that  the 
hardware  and  software  implementations  used  are  perform- 
ing satisfactorily. 

2.  The  log-normal  distribution  seems  to  be  the  number  one 
candidate  for  most  acceptable  fitting  distribution  for 
data  whose  amplitudes  are  proportional  to  the  square 
root  of  intensity.  There  seems  to  be  some  promise  that 
an  analytic  distribution,  such  as  the  log-normal  or  the 
gamma,  can  be  contaminated  (see  section  3.2.3)  so  as 
to  produce  a distribution  whose  fit  is  acceptable  to  a 
statistical  test  such  as  a chi-square. 

3.  The  clutter  data  for  pack  ice,  bare  mountain,  city  and 
field  provide  a values  which  are  in  reasonable  agree- 
ment with  the  small  amount  of  L-band  clutter  data 
presently  existing  in  the  literature. 

4.  When  the  resolution  cell  size  is  increased,  i.e.,  when 
the  data  are  averaged  over  range-azimuth  rectangles  of 
increasing  size,  the  standard  deviations  of  the  result- 
ing distributions  of  averaged  data  decrease  correspond- 
ingly. The  rates  of  standard  deviation  decrease  vary 
with  incidence  angle  as  well  as  with  the  type  of  terrain. 

The  following  recommendations  are  made  with  regard  to  additional 
work  in  clutter  statistics  by  ERIM,  RADC , or  other  organizations. 

There  is  now  a rather  complete  computer  program  library  to  generate 
clutter  statistics  which  will  make  the  analysis  of  clutter  data  more 
efficient  in  the  future. 
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1.  Since  some  questions  remain  on  the  behavior  of  the  clutter 
data  at  the  smaller  depression  angles,  the  attempted  cali- 
bration at  these  angles  should  be  repeated  using  the  assump- 
tion thac  the  RCS  of  the  corner  reflectors  imaged  is  correct. 
This  technique  avoids  the  problem  of  knowing  the  antenna 
pattern  and  it  does  not  restrict  calibration  to  depres- 
sion angles  where  there  are  corner  reflectors,  i.e., 

26°,  27°,  30°,  42°,  and  44°. 

2.  The  radar  simulation  derived  by  DSA  and  described  in 
Section  4 of  this  report  was  used  to  provide  magnetic 
tape  data  for  use  in  testing  signal  processing  techniques. 
However,  with  the  more  extensive  data  now  available,  more 
realistic  radar  clutter  data  simulation  could  now  be  per- 
formed. It  is  recommended  that  work  continue  in  the  areas 
of  appropriate  statistics,  generation  of  clutter  maps 
and  on  determination  of  the  spatial  correlations  of  actual 
data. 

3.  In  order  to  analyze  the  data  at  coarser  resolutions  cor- 
responding to  specific  radar  designs,  long  (160  km)  strips 
of  data  can  be  digitzed  and  analyzed.  The  160  km  strip 
could  be  digitized  with  100  m range  resolution  and  a variety 
of  cross  range  resolutions.  This  would  lead  to  an  investi- 
gation of  clutter  statistics  as  a function  of  terrain  types 
along  a continuous  160  km  strip.  The  statistics  thus  derived 
would  apply  directly  to  the  design  parameters  of  the  space 
based  radars  being  considered.  The  software  developed 

on  the  present  program  can  be  used  to  derive  the  required 
statistics. 

4.  There  are  a number  of  tasks  which  would  improve  accuracy  and 
establish  additional  credibility  to  the  clutter  statistics. 
These  would  concentrate  on  analyzing  and  interpreting  the 
results.  The  present  calibration  technique  can  be 
analyzed  to  establish  achievable  accuracy.  The  computer 

-3 

programs  can  be  extended  to  normalize  out  the  R term 
and  the  antenna  elevation  gain  pattern  from  the  digitized 
data.  This  presently  restricts  us  to  looking  at  statis- 
tics along  track  and  over  small  elevation  angles. 

Other  statistical  distributions  can  be  examined  and 
compared  including: 

a.  contaminated  log-normal 

b.  contaminated  gamma  (or  chi-square) 

c.  Beckmann  Distribution, 

d.  Rice  Distribution  (or  chi),  and 

e.  Hoyt  Distribution. 
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In  addition,  it  may  also  be  useful  to  attempt  a nonparametric 
estimation  of  a density  function.  Also,  other  statistical 
tests  can  be  studied  to  further  define  a "good  fit".  Other 
tests  besides  the  chi-square  which  can  be  examined  are: 

a.  Komolgorof f-Smirov  Test, 

b.  Snedecor  Test, 

c.  Pains  Test. 

The  spatial  correlation  properties  of  the  data,  as  well 
as  spatially  varying  properties  of  the  distributions, 
can  also  be  examined . 

5.  In  order  to  derive  a parametric  model  for  clutter,  one  can 
study  the  relationship  of  radar  cross-section  per  unit  area 
( a0 ) as  a function  of  wavelength,  polarization,  aspect 
angle,  terrain  type,  etc.  A model  then  could  be  used  in 
analytic  or  simulation  radar  performance  evaluations. 

6.  The  present  work  could  be  extended  to  X-band  (using 
available  data)  and  to  sea  clutter  at  X-band  and  L-band. 

7.  Cross-polarized  data  can  be  used  to  investigate  whether  the 
use  of  a cross-polarized  channel  can  result  in  a significant 
decrease  in  clutter  return.  This  would  include  a review 

of  limited  data  in  the  L-  and  X-band  cross  polarized 
channels  of  the  ERIM  radar  and  comparison  to  the  same 
polarization  channels. 

8.  If  future  programs  justify  extensive  additional  data 
collection,  then  we  would  recommend  that  adequate  amplitude 
calibration  instrumentation  and  procedures  be  added  to  the 
X-L  radar  to  ensure  accurate  calibration. 

9.  The  ERIM  X-  and  L-band  dual  polarization  radar  used  to 
collect  the  imagery  used  in  this  program  is  in  the  process 
of  being  transferred  into  another  aircraft.  The  new 
aircraft  (a  convair  580)  will  be  used  in  a pilot  program 
to  determine  microwave  sensor  requirements  for  Canada. 

This  program  will  generate  a minimum  of  15,000  line  miles 
(54,000  square  miles)  of  radar  imagery.  Included  in  the 
data  set  will  be  agricultural  areas,  wildlife  preserve 
areas,  and  forested  areas.  Also  included  will  be  ice 
studies  from  both  coast  lines  as  well  as  the  artic 
regions.  The  open  ocean  and  coastal  areas  will  also  be 
imaged . 

The  imagery  collected  during  this  program  (July  78  through 
April  79)  will  be  available  to  ERIM  for  clutter  statistics 
analysis . 

There  is  at  the  present  time  no  funding  for  adding  amplitude 
calibration  to  the  X-L  radar.  The  value  of  the  data  collected  could 
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be  enhanced  by: 


1.  Adding  amplitude  calibration 

2.  Conducting  an  error  analysis  on  the  present  system  to  put 
error  bounds  on  the  data. 

This  system  will  remain  available  to  ERIK  and  its  sponsors  for  the 
duration  of  the  program  through  the  normal  channels. 
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